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' Strongly correlated two dimensional electrons are believed to form a spin liquid in some regimes of 

density and temperature. As the density is varied, one expects a transition from this spin liquid state 
' to a spin density wave antiferromagnetic state. In this paper we show that it is self-consistent to 

assume that this transition is second order and, on this assumption, determine the critical behavior 
of the 2pF susceptibility, the NMR rates Ti and T2 and the uniform susceptibility. We compare our 
results to data on high Tc materials. 



I. INTRODUCTION 



I High temperature superconductors may be created adding carriers to magnetic insulators. At low dopings the 
compounds have long range magnetic order at T = 0; at high doping they do not. Therefore, a T — magnetic- 
non-magnetic transition must occur when the carrier density exceeds a critical value. The idea that some of the 
. anomalous properties of high Tc superconductors are due to their proximity to this quantum phase transition has 
' attracted substantial recent interest. 

The properties of the transition depend on the ordering wavevector and on the nature of the disordered phase. Sev- 
'■^ ' eral different possibilities have been studied in some detail including antiferromaftoet-singlet transitions in insulating 
^ ' magnetsEI and ferromagnetic and antiferromagnetic transitions in Fermi liquidsacl. Here we consider an important 
O • case which has not so far been discussed in the literature, namely, that the disordered phase is a "spin liquid" and 
the ordering occurs at the wavevector |Q| — 2pp. 

By "spin liquid" we mean a liquid of charge zero spin 1/2 fermion excitations coupled by a singular gauge field 



intecaction; in the ground state the fermions fill a large Fermi sea which occupies a substantial part of the Brilloui 
zonaj. The spin liquid model has been argued to describe the normal phase of high temperature superconductorsf 
' By "2pF^^ we mean a wavevector which connects two points on the Fermi surface with parallel tangents (See Fig. l)t 
For a circular Fermi surface any vector Q of magnitude 2pF connects two such points. One important motivation for 
studying the 2pp case is the high Tc superconducting material La2-xSrxCuOi, in which strong magnetic fluctuations 
have been observedEl; the fluctuations are peaked at an x-dependent wavevector Q(a;) which is clairpjed to be a "2^^" 
wavevector of the Fermi surface calculated by standard band-structure techniques for this materiafl. 

The 2pp spin density wave transition in a spin liquid is different from the 2pp transition in a conventional Fermi 
liquid because the gauge field interaction leads to divergences in the fermion response functions at wavevectors |Q| 
near 2pp. In this paper we calculate the exponents and the scaling functions characterizing the 2pp transition in the 
spin liquid. Our starting point is a Hamiltonian, H = HpL + Hgauge, where Hp^ describes fermions moving in a 
lattice and interacting with each other via a short range four fermion interaction W: 

HfL =^eip)cl^aCp^a + W ^ c], „Cp+q,a4,_^Cp'_,,/3 (1-1) 
pa p.p' ,q,a,(} 



and Hgauge describes the gauge field and its coupling to the fermions; it is discussed in detail in the literatureaa and 
below in Section H. 
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Figl 



FIG. 1. Sketch of Fermi surface and important wavevectors. The Fermi surface shown here is similar to that claimed to be 
appropriate to Lai.^Sr .liCuOi ■ The vector Q connects two points on the Fermi surface. It is assumed that the tangent to 
the Fermi surface at one end of the vector Q is parallel to the tangent to the Fermi surface at the other end. We parametrize 
the vector k by the angle 6 shown on the figure and the magnitude fey shown on the figure as the solid line part of k. 

It has recently been shown that two possibilities arise, depending on the strength of the fermion-gauge-field 
interactioro. If this interaction is weak, a critical value, Wc, of the short-range four-fermion interaction separates a 
disordered phase with non-divergent spin fluctuations at < Wc from an ordered phase aXW > Wc- In this paper we 
evaluate the exponents characterizing the transition at W — Wc in the weak gauge-field coupling case and show how 
they depend on the strength of the fermion-gauge-field interaction. If, however, the fermion gauge-field interaction 
is sufficiently strong, theii,the 2pF spin susceptibility diverges as T — )■ for arbitrary W < Wc, although there is no 
long range order at T = 03. We have not studied the transition at W = Wc in the strong gauge-field interaction case, 
but we give the exponents characterizing the W < Wc phase and discuss the physical consequences. 

We assume that transition at W = Wc leads to a spin density wave state with long range order with wave vector D 
and that this transition is of a second order- If |Q| 7^ 2pi?, the gauge field does not modify the fermion susceptibilityD, 
so we expect previously developed thcoryclu to apply. We now outline our approach to the |Q| ~ 2pp case. Because 
we expect the physics in this region to be determined by the exchange of spin density fluctuations we use a Hubbard- 



Stratonovich transformation to recast Eq. (1.1) as a theory of fermions coupled to spin fluctuations Sqi 



H' = ^e{p)cl,aCp.a+ g ^ cl^^aaf3Cp+k,l3S-k (1-2) 
pa p,k,a,l3 

+Hgauge + ^ S| 

k 

Here g is a fermion-spin fluctuation coupling constant derived from W; in weak coupling = W. One may study 
H' as it stands or one may integrate out the fermions completely, obtaining a theory of interacting spin fluctuations 
which is described by the action 

A[S] ^Y.^^''^^oHk,^)S-u..k (1.3) 

uJi^ki i 

Here Suj,k represents a spin fluctuation of Matsubara frequency uj and wavevector k. xo is the susceptibility and U is 
a four-spin-fluctuation interaction proportional to g*. We shall derive and interpret this action in more detail below. 
This action is difficult to treat even for a Fermi liquid without the gauge fleld interaction because xo ^md U diverge 



2 



as T, u; — > jQi-aind k Q. The gauge field causes additional singular renormalizations of xo, U, g and the fermion 
propagatoroE^. 

In this paper we present and justify a a self-consistent one-loop approximation method for extracting physical results 



from the formally divergent theory. We supplement this treatment by an analysis based on Eq. (1.2) of the effect of 
the spin fluctuations on the fermions. We find that the susceptibility is less singular than the susceptibility obtained 
for transitions with Q ^ 2pp. This is because the 2pp singularity of the fermions leads to a long ranged (RKKY) 
interactions; this weakens the singularities associated with the transition for the same-reason that long ranged dipolar 
interactions lower the critical dimension of classical models of ferroelectric transitionO. 

The self-consistent one loop approximation succeeds in the spin liquid case because the fermion-gauge field interac- 
tion makes the singular part of the fermion response symmetric in the variable |k| — 2pp. This approach however fails 
for the conventional Fermi liquid because in this case the singularities in the fermion polarizability are not symmetric 
in the variable |k| — 2pp as discussed below. We have found that this asymmetry implies that fluctuation effects do 
not permit a second order transition at |Q| = 2pp in a Fermi liquid. We will present a detailed treatment elsewherdiS. 

The outline of this paper is as follows. Section II reviews the relevant theory of the spin liquid. Section III is devoted 
to the transition in the spin liquid for weak fermion gauge-field interaction. Section IV discusses the properties of the 
small- critical phase occurring when the gauge interaction is strong enough that the 2pp susceptibility diverges as 
T ^ 0. Readers uninterested in the details of the derivations may proceed directly to Section V which contains a 
summary of the results and a conclusion. 

II. GAUGE THEORY OF SPIN LIQUID 

In this section we review the properties of the "spin liquid" regime of the t-J model, which has been argued to 
describe the normal state of the high-Tc superconductorsQ. In the spin liquid regime the elementary excitations are 
charge e, spin "holons", charge 0, spin 1/2 fermionic "spinous", and a gauge field which mediates a long ranged, 
singular interaction. For temperatures above the Bose condensation temperature of the holons, the holons have a 
negligible effect on the magnetic properties, which are determined by the properties of the degenerate Fermi gas of 
spinous (which we call a spin-liquid), coupled by the gauge interaction and by a short-range four- fermion interaction 
W. 

The properties of the spin liquid have been studied by many authors@@'i. The fermion propagator has been found 
to have the form 

G=^y5 (2-1) 

iuj^' ^2/3 _ +p\/pq) 

instead of the Fermi liquid form G — — w(p|[ +p'\^/pq)]~^ . Here vp is the Fermi velocity, po is the radius of curvature 
of the Fermi line and p\\{p±) are momentum components normal (tangential) to the Fermi lin£ as measured from the 
points ±Q/2 and ujq is an energy scale of order ppvp which is defined more precisely in Ref.u. 
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Fig 2 



FIG. 2. (a) Diagram yielding the nonanalytic momentum and frequency dependence of the susceptibihty xo- The solid 
lines are the fermion propagators. The shaded triangles are vertices which have power law divergences, (b) Diagram yielding 
the nonanalytic momentum and frequency dependence of the four spin fluctuation interaction U. 



Recently, we have shownH that the fermion-gauge-field interaction causes a power law divergence in the vertex r2pf 
which is shown as a black triangle in Fig. 2a. Specifically we found, 



2pf 



We were able to calculate a in the limits N 
found 



a 
a 



max 



(—) (- 

\uja J ' \ u. 



3/2 



(2.2) 



— > oo and — > 0, where TV is the spin degeneracy of the fermions. We 
l/2iV N^oo 



Extrapolation of these results to the physical case N — 2 gives the estimates 1/4 ^ tr ^ 3/4. The combination of the 
change in form of the Green function and the divergence of T2pp has profound effects on the polarization operator 
TUlo, k) (which is shown in Fig. 2a) and in the four spin fluctuation interaction U (which is shown in Fig. 2b, cf. Eq. 

&)■ n 

If 1/6 < cr < 1/3, we foundEl 



n(cc;,fci|) = no- 



H 

UJQ 



CkRe 



l-3a' 



(2.3) 



Here and are constants of order unity. For a < 1/6, the term involving fc|| has the same sign as /c||. Thus 
the maximum of n|-|Smd therefore the ordering wavevector must occur at a |Q| different from 2pp and we expect 
previous treatmentaalj of the critical behavior to apply. For cr > 1/6 the fc|| term is always positive and a second order 
transition at |Q| ~ 2pp is possible. We restrict our attention to cr > 1/6 in the rest of the paper. In this case the 
nonanalyticity of i?efcj|^'^'^ affects only the magnitude of c^, that is not important for the scaling arguments which we 
will present. We henceforth write the fcii term as 



Ck 



k\\VF 

Wo 



1-3(7 



On the other hand, if ct > 1/3, 11 diverges as k« 0. We founc 
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The full susceptibility x is obtained by combining the irreducible bubble 11 with the short-range four fermion vertex 
W . We have showntl that the gauge-field interaction renormalizes a sufficiently weak initial Wi to zero. Therefore we 
expect a T = transition as W is varied through a critical value Wc- In the weak coupling case cr < 1/3 the fact that 
n(a; ~ 0, fc|| ~ 0) is non-divergent implies that the usual RPA formula 

is the correct starting point of a theory of the transition. We discuss the transition occurring when W = Wc = 
n(0, 0)""'^ in more detail below in the next section. 

In the strong coupling, a > 1/3, case we still expect a transition when W exceeds a critical value. However 



the divergence of 11(0,0) implies that the RPA formula (2^) is not correct. For W < Wc, x(a;,fc) = n(w, /c) with 



corrections of order the product of W and a positive power of frequency or fcy . The T = critical point separates a 
W < Wc phase which has power-law spin correlations from a.W> Wc phase which has long-range order. 

III. WEAK FERMION-GAUGE-FIELD COUPLING 

There are two possible approaches to the spin density wave transition a.t W — Wc- One is to integrate out the 
fermions, obtaining an effective action, .4[S'], describing the spin fluctuations, S^j^k, and then to analyze this action. 
If Q ^ 2pi? all terms in the effective action arc finite in both .the Fermi liquid and in the spin liquid, and the action 
may be treated by standard renormalization group tcchniquesa. As Q 2pp, however, the coefficients of the quartic 
and higher order terms diverge, implying that standard KG techniques cannot be used and that the effect of higher 
order terms must be investigated. Another approach is to investigate the effect of the spin fluctuations on the fermion 
propagator. Here we consider both approaches. 

Taking the first approach we proceed in four steps. First, we integrate out the fermions obtaining the action A[S] 



given in Eq. (F3). Second, we truncate the action, retaining only the quadratic and quartic terms. Third, we solve 
the truncated action in the self-cpiisistent one loop approximation. This approximation has been extensively used to 
study three dimensional magnetgll3 and gives the same results as the renormalization group above the upper critical 
dimensionEl. Fourth, we show using the formalism of the effective action A[S] that corrections to the self-consistent 
one loop approximation lead at most to logarithmic corrections to the self consistent one loop results. Fifth, we 
confirm this result using the formalism of the spin fluctuations interacting with fermions. Finally, we discuss the 
physical consequences of the theory. 

A. Derivation of truncated action 



We begin by evaluating the coefRcients xo and U in action ( |l.3| ). We shall be interested in momenta close to the 
momentum Q at which xo is maximal. For wavevectors near Q the momentum and frequency dependence of x and 
U are non-analytic and controlled by Fermi surface singularities. It will be convenient to parametrize the momentum 
k in terms of a magnitude and an angle 9 as shown in Fig. 1. For each angle 9 there is a momentum Pf{9) such 
that 2pp{9) spans the Fermi surface. Note 2pp{9 = 0) = Q. We define fcy to be the difference |k| — 2pp{9). These 
definitions are generalizations for the noncircular Fermi surface of coordinates which are the convenient choice for a 
circular Fermi surface. 

The fermion contribution n(u;,fc||) to the inverse susceptibility Xo^ can be calculated by summing all diagrams 
which are irreducible with respect to the fermion-fermion interaction and have two external S^j^k legs. This sum has 
contributions from short length scale processes which give n(ti;,fc||) an analytic dependence on fc|| and 9 and also 
contributions from Fermi surface singularities, which lead to a nonanalytic dependence of H(a;, fc||) on fc|| and uj. The 



Fermi surface singularities come from the diagram shown in Fig. 2a which leads to the expression (2.3). Thus we get 
the susceptibility 



(3.1) 



where g is the fermion-spin fluctuation couphng. The coefficients Cui and Ck are of the order of unity, they are sensitive 
to the details of the band structure and the momentum dependence of the interaction. The coefficient Kq is determined 
by difference of the interaction from its critical value. Note that we are using Matsubara frequencies so that x is 
purely real. 

The most singular contribution, Using, to U is given by the diagram shown in Fig. 2b. It diverges if the reduced 
momenta and frequencies of all four legs are zero. For our subsequent calculations we shall need to estimate 

the asymptotic behavior when the momenta and frequencies on the external legs satisfy 



•|| 



W{, > CJc ~ 



and 9a ~ Oh', 0c ~ Od but with 9a — 9c arbitrary. Here a,b,c,d are any permutation of the legs shown in Fig 2b. By 
evaluating the diagram in Fig. 2b we find that in this limit the vertex is 



Um,u;;9i,92)^^ 



1 



/H\4'^+2/3 f vFP%ie^-02V \ 



60-+1 



6(T + 1 



(3.2) 



The equation (3.2) is not strictly correct at arbitrary momentum and frequency but gives the correct asymptotics of 



the vertex, which are all we need. The result may be motivated by the following arguments: the scaling loq' ' lo 
VFk\\ ^ vpippOY /po comes from the structure of the expanded fermion Green function ( |2.lD ; we believe that this 
is the correct generic scaling in this problem. In the case of the circular Fermi surface the dependence on the large 
angular separation (9a — 9c) follows from rotational invariance; a noncircular Fermi surface can still be mapped locally 
to a circle, proving this result in the general case. In writing these equations we have included the renormalization 
of ^2pF due to the fermion-gauge-field interaction. The renormalization of the fermion vertex due to the gauge 
field does not enter our discussion of the critical phenomena. To see this, note that the renormalization is due to the 
logarithmic divergence of the diagram shown in Fig. 3. In the present case the interact ion is sharply peaked, so the 
logarithm is confined to momenta and frequences less than the scale set by k (see Eq. (3.1)), and therefore does not 
affect our tre at men t of the critical point. The problem is thus to investigate the self-consistency of the theory defined 
by Eqs. (O), dp) and (O). 



FIG. 3. Diagram yielding renormalization of the spin-fluctuation contriubtion to the four fermion vertex W by gauge field. 
Here the solid lines are fermions, the dashed line represents the gauge field, and the wavy line represents the spin. 



B. Self-consistent one-loop approximation 



We treat the theory defined by Eqs. (1.3,3.1,3.2) in the self-consistent one loop approximation shown diagrammat- 
ically in Fig. 4, i.e. we require that the full susceptibility obeys the equation 
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X{k, uj)-' = Xo"'(fc, ^) + / E 



(3.3) 



where xo is given by (|3.l| ). Using (3.1) and ( p.2| ) we see that the frequency or momentum derivative of the integral 
(|3.3[) is infra-red divergent. Estimating this integral we get: 



-in ^ -In ^ 4 l ^oPo {HV^^ , 
X {k,u;)^Xo ik,u.)-9 [-) + 



VF\k\\ 
Wo 



(3.4) 



which has the same order of magnitude and the same frequency and momentum dependence as the bare Xo ^ • Here as 
in Eq. (3.2) the formula gives only the correct asymptotics when one of the variables (wg^'^a;^/"^, u_fA:|| , vf{pfQY /po, 
vpK^ /po) is much larger than the others. 
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Fig 4 



FIG. 4. Diagrams defining the self-consistent one-loop approximation for the susceptibility x (thick wavy line) in terms of 
the interaction U (shaded box) and the bare susceptibility xo (thin wavy line). 

Due to the four spin interaction U the cutoff, k']^(T), is renormalized as compared to Kq. At a finite temperature 
the sum over rj in (|]^) is taken over Matsubara frequencies rjn — 2'KnT . The temperature dependent part of the 
correction is dominated by the first term in this sum. Estimating it at zero external momenta and frequency we get 
the leading contribution to 



Pi 



pI 



(3.5) 



Note that x~^{S) ~ 6*^; the singular interaction does not change the power law in the angle dependence. This may 
be most easily seen via a reductio ad absurdum. Suppose the fou r-spin-wave interaction had led to an exponent less 
than two for 6. Then th e an gular dependence of xo in Eq. (3.2) could have been neglected. However for a circular 
Fermi surface Eqs. ( p.2| ,3.3) are rotationally invariant (apart from terms due to Xo)j ^-nd can therefore lead to no 
angular dependence at all, in contradiction to the hypothesis of angular dependence with an anomalous exponent. 
The same argument applies to a nonspherical Fermi surface, because it may be mapped into a spherical one, with 
errors of order 0'^ . We conclude that in general the self-consistent one- loop equations cannot produce a ^-dependence 
different from 9^. 



C. Marginality of Higher Order Corrections 



We now argue that corrections t o the se lf-consistent solution (3.1) do not change the momentum, frequency and 
temperature dependences of Eqs. (3J,^). There are two kinds of corrections: those arising from higher orders of 
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perturbation theory in the S"* couphng u sing the truncated action, Eq. (2.1), and those arising from higher order 
S^,S^... nonlinearities omitted from Eq. (L3). We first consider the renormahzation of the vertex U at the second 
order in U. The corresponding diagram is shown in Fig. 5. This diagram is infrared divergent and may therefore 
by evaluated by writing the integrals with zero external momenta and frequencies and then cutting off the resulting 
divergence by the largest of the external momenta or frequency. This gives 



du{n,q,ei,92) 



g^PF 



dujdOdkw 



aj>0 fc|| >q 



, s 2/3-2. 



l-3cr 



4(7+2/3 



, 6cr+l 



1 












■(s) 









4(7+2/3 



+ 



6(7+1 



g_ l^oPa 



1 



vf 



6(7+1 



, 6(7+1 



(3.6) 



This estimate shows that higher order corrections in the 5^ interaction do not change the power law dependence of 
the mode coupling interaction, and therefore cannot change the powers coming from the solution of the self-consistent 
equation. This argument of course does not rule out logarithmic corrections to ■ 









+ 












FIG. 5. Diagram giving leading renormalization of full four spin interaction U (shaded box) in terms of bare U (open box) 
and spin fluctuations (wavy lines). 



The effect of the higher order nonlinearities {S^, S^, etc) may be analyzed similarly. Let us compare the contribution 
to the susceptibility from one bare 5^"+^ vertex to the contribution from one bare 5^" vertex. The former contains 
one more integral over (a;,fc||,0) and one more factor of the susceptibility x- The diagram for the bare vertex 
S'^""'"^ contains also two more fermion Green functions. Performing the extra integration we see that the additional 
singularity coming from the extra fermion Green functions in the bare vertex and extra factor of susceptibility is 
precisely cancelled by the phase volume. 

These results imply that the theory is marginal, in the sense that higher order interactions give the same infrared 
behaviour as lower order interactions. The marginality may also be demonstrated by a scaling argument. If we rescale 
momenta and frequencies via lj' = bu>, fcy = b^/'^fcy, 9' — b^^^9, then we must scale the field S'^ — b^^/'^'^'^Sq to keep 
the quadratic term in the action invariant. The scaling dimension of the Sf- vertex is then 6(-8/3+2(7)n (^^^^^ ^^le 
fields) times fe-4n/3+2~2n(7 (|fj.Qm \}^Q vertex) times 62(2n-i) (^fj-Qni the integrals). Adding the powers we see that the 
total scaling dimension of the vertex is zero, so all interactions are marginal in the renormalization group sense. 
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Thus, we conclude that higher order effects of the spin wave interaction do not chan ge th e exponents characterizing 
the divergence of the susceptibility or the interaction between spin fluctuations (3.2). 



D. Effect of spin fluctuations on fermions 

We now c onsi der the alternative approach of estimating the feedback of the interaction mediated by the spin 
fluctuations ( |3.lD on the fermions. The lowest order contribution to the fermion self-energy S is 

(3.7) 



S(e,p) - .g2 / G{e',p')x{e-e',p-p')[r^p+p,y,^^,^,,y,{e-e\p-p')]'{dp'de') 



We are interested in the leading frequency dependence of E; the main contribution to this comes when p' is on the 
Fermi line. We then find the angle at which x is peaked; this turns out to be the point which is located symmetrically 
opposite to p: \9' — 6 + tt\ ^ (k^ + 46*^) 2(i-3<t) where 6' and 9 are the polar coordinates of the points p' and p. 
Estimating the contribution of this region we get: 



S(e,e) = 



PF ' 



(3.8) 



where ce ~ 1. Over most of the Fermi line, this self-energy is small compared to the fermion self-energy ujy^e^/^ due 
to the gauge field interaction. The spin fluctuation contribution becomes important only in a small region 16*1 < 6* 
near the points connected by the wave vector Q, with 9* given by 



(l-3cr) 

= (e/t^o) ^ 



At a finite temperature T the typical energy of the fermion is T, and it should replace e in (3.9) 



r{T) - (T/^o)^ 



(3.9) 



(3.10) 



remains of the order oi uj]/^ e^/^ . 



Since 9*{T) > n/pp (cf. Eq. (gj)), in the region 9 < 9*{T) the self-energy (E 

We do not have a reliable expression for the self-energy S](e, 9) in the region \9\ < \9* \ at T. We can, however, 
argue that fermions in this region make only a marginal contribution to physical quantities such as polarization 
operator n(w, k). The analytic expression corresponding to Fig. 2a for 11 is 



n(a;,fc||,0) = J dpdeT%^G{p + k,e + u;)G{p,e) 



(3.11) 



Here k is a vector parametrized as shown in Fig. 1 by variables As a function of angular variable 9, Il(u;,9) 

has a smooth maximum of the width 9'{lj) ^ {to IujqY/'^^'' ^ 0*(t j). Thus all physical quantities are controlled by 
spin fluctuations with 9 > &'{^)- Using Eq. (2.1) for G and ( |2.2| ) for T2pp we see that the integral over momenta 
is dominated by fermions with momenta near the Fermi surface which have angles very close to 9 (the dominant 
contribution comes from the range of fermion ang les 9': \9' - 9\ ^ {uj/uaf/^ < 9). Thus for k\\ and 9 which 
are important for the susceptibility, for the self-energy S(e) or for the T-dcpendence of k, the polarization bubble is 
controlled by the fermions aX 9 ^ 9* . Thus, we believe that the fermion-spin-wave interaction leads only to corrections 
of the order of unity to physical quantities. Certainly, these arguments are based on power counting and can easily 
miss logarithmically large contributions; however, we have verified that the leading vertex and self energy corrections 
do not contain logarithms. 



E. Physical consequences 



We now discuss the physical content of the results. We first note that the fermion-gauge field interaction has 
two effects on the polarizability near Q — 2pF'- it changes the form of the nonanalyticity at w = and Q = 2pp 
(introducing the exponent a) and it washes out the nonanalyticity associated with the lower boundary of the particle- 
hole continuum at w = 2vpk\\. Therefore the scaling form, eq. (3.1), gives the correct resu lt fo r the imaginary part of 
the susceptibility, x", which may be measured in neutron scattering. Using this and eq. (2.5) gives 
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(l-g2n(5,a;))2 
The NMR Ti relaxation rate is therefore 

/mn(/c|| , uj) 



X[q,'^) = T, „?TT^.. , .^^2 (3-12) 



-i- = A'g^ [ d^pFde lim ^"""^^"'""^ [x(0, , 0)]^ (3.13) 

Here A is a constant proportional to the hyperfine couphng. We can neglect the weak dependence of 11 on 9 since the 
singular dependence on 6 comes only via x- For T > and fc|| small (fc|| < pp{T/ J)"^^^) we have 

lim /""(^-"="'"^ = / G(.,p + Q/2)G(.,p-g/2)[r(«)(0,Q)]^ - 



0)^0 



2Tcosh2(e/2T) 



CTj^i^Y (3.14) 



while for larger fcii it decreases as 



/y2/3 1/3^^/2 

lim /mn(fc||,w) - /mn(0,w) -2— (3.15) 

uj-^a y VF\k\\\ J 

The divergence of /r7in(fc|| ~ 0,lu)/uj is the usual Kohn anomaly modified by the fermion-gauge-field interaction. 
This interaction has two effects. First, the increased fermion damping proportional to e^/^ weakens the T dependence 
from T~^/^ to T~^/^ but does not change the q-dependence. Second, the extra vertex correction then strengthens the 
T dependence to y-(2o-+i/3) ^^^^ ^]-^g g-dependence to \q — 2pi?|~^'^'^+^/^^ . If cr > 1/6 as we assume throughout the 
strengthened g-dependence leads to a divergence even far from the critical point. 

If the system is tuned to the critical point, then the T dependence of x becomes important and we obtain 

(3.16) 

If (7 < 1/3, the critical contribution is more singular than the background contribution. However as a is increased to 
1/3, the difference between the critical contribution and the background contribution disappears. 

One may calculate the T2pcate similarly. The electronic contribution to the NMR rate is related to the real 
part, x', of the susceptibilitytJ: 




where a is the lattice constant. If o- < 1/9, the critical contribution to ^ is 

2 ^ ag^y pAt) 

Since we have assumed a > 1/9, is not divergent as T — > 0. 

Proxinwtjy to the antiferromagnetic transition has also an effect on the uniform susceptibility. We have shown 
elsewherellJ that the leading low-T behavior of x(q, 0) in the limit q ^ is given by 



Xu 



lim x(q,0) - Vi?(fc,;/)x(fc,i^) (3.17) 



where the coefficient D is given by the diagram in Fig. 2b with two of the spin fluctuation vertices replaced by vertices 
coupling the fermions to the external magnetic field. We denote this coupling to the external magnetic field by gg. 
In particular, terms arising from triangular vertices coupling two spin fluctuations to one small-g external field lead 
only to terms proportional to integer powers of T, w hich are much smaller than the terms we keep. The calculation 
of D is very similar to the calculation of U (Eq. ( |3.2| )) except that the coupling D has two small q vertices and only 
two large Q vertices r2pp ; the degree of singularity therefore differs from that of U and the dependence on 9i — 62 is 
absent. We find 
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1/2 

Po 



Substituting this expression for D and Eq. {"iA) for x into Eq. (2.13) we find at g 



3cr+l 



— that 



Xu ^ const + Dq (T/ cjq) 3 " 



(3.18) 



(3.19) 



where £>o is a constant of order glpp/ig^vp)', Do is positive if n/pp > {T/uJaY^'^^'' ^^'^ is negative if n/pp < 
{T /ujqY/^^" . The rel ation between the sign of Do and the magnitude of k conies because there are two sources of 
T-dependence in Eq. ( [3.17| ): one is the T-dependence of the cut-off k; the other is the discreteness of the frequency 
variable. 

If k{T = 0) > then a crossover occurs. For T^-" > K, the result (|3.19|) holds. For lower T , the renormalization 
of the 2pp interaction by the gauge field shown in Fig. 6 becomes important and the T dependence becomes weaker. 
We discuss the T-dependence in the g < gc regime in the next section because it is important for cr > 1/3. 




Fig 6 



FIG. 6. Diagram yielding leading T-dependence of the uniform susceptibility for the case of strong gauge-field-fermion 
coupling. The solid lines are fermion propagators and the heavy dot is the short-range four fermion interaction dressed by the 
gauge field. 



IV. STRONG GAUGE FIELD INTERACTION a > 1/3 



If the exponent a arising from the fcrmion-gauge-field interaction is greater than the critical value 1/3, then the 
fermion spin susceptibility diverges at q — 2pp and lu = for all g < gc, leading to a phase with power-law spin 
correlations at T = 0. At g = gc & transition will occur to a phase with long range order. We have not yet formulated 
a theory of this transition. In this section we outline the physical consequences expected for the g < gc phase. 

Fo r a > 1/3 the Kohn singularity at q = 2pp becomes nonintegrable, so x' and x" /'^ diverge at T = as shown in 
Eq. (2.4). The power law describing the divergence is ^-independent. The NMR relaxation rates are therefore given 
by summing the appropriate combinations of susceptibilities over ku . We find 



1 



A 



1/2 1/2 / rri 

2PfPq' ujq' I T 



5/2 



If a < 1/2, the rate ^ is non divergent and if cr > 1/2, 



1 



UJQ 



-2a 



1-2(T 



(4.1) 



(4.2) 
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Here a is the lattice constant. We now consider the uniform susceptibiHty x{q ~^ 0,0). It has contributions from 
diagrams involvingpthe four fermion vertex W. The leading diagram is shown in Fig. 6. Because, as we have shown 
in a previous papercl, W is renormalized by the gauge-field, xu acquires an anomalous temperature dependence. The 
renormalization of W is cut off by the largest of the temperature, the energy of any fermion line, the scaled momentum 
(wpil)^/^ of any ferrpion line, or the difference from tt of the angle 9i — 62 between the momentum of the incoming 
particles, leading too 



W ^ Wo max 



3/2 

, (01 - 02 - 7r)3 



(4.3) 



where the exponent (3 is iV-dependent and is not simply related to cr. As iV ^ 00, I3{N) — ^ — j^; as N ^ Q 
j3{N) ~ ci^/N with ci > 0. 

To calculate the temperature dependence of x(9 — > 0, 0) we insert the vertex given in eq. ( [4.3| ) into Fig. 6, and then 
determine the phase volume in the (e,p), ,p') integrals in which the T-dependence of W is important. We find 

Xu= const + D'^{T/uof+f' (4.4) 
where Dq is a constant of order ^^^j^^ W, whose sign is positive for repulsive W and negative for attractive W . 



V. CONCLUSION 



We have determined the scaling behavior of the spin susceptibility near an antiferromagnetic critical point at a 
"2^;^" wavevector of a spin liquid. By a "2^;^" wavevector we mean one which connects two points on the Fermi 
surface with parallel tangents. An example is shown in Fig. 1. We distinguish between very weak, weak and strong 
fermion-gauge-field coupling. The cases are defined by the value of the exponent a appearing in Eq. ( |2.2| ). This 
exponent depends only on the fermion spin degeneracy, N] we do not know how to calculate it analytically for the 
physically relevant case N — 2. For estimates of a, see ref.O. In the very weak coupling case, a < 1/6, a second 
order transition at |Q| — 2pp is impossible because the maximum of x is at |Q| 7^ 2pi?. In the weak coupling case 
l/6<(T<l/3 the critical point occurs at T = when a short-ranged interaction W equals a critical value Wc- In 
the strong coupling case cr > 1/3 there is a critical phase with power law correlations at T = for < W < Wc and 
an additional transition, which we did not study, a.t W = Wc- 

In weak and strong coupling cases the 2pp singularities of the fermions lead to physically important effects. First, 
the scaling is anisotropic. The dependence of the susceptibility on wavevectors fc|| which are parallel to the ordering 
wavevector Q involves a different exponent than does the dependence on wavevectors k± perpendicular to Q. The 
fcx e xpo nent takes the conventional Ornstein-Zernike value 2 while the fc|| exponent is always between and 1 (cf. 
eq. ( p.l| )). The small value of the fc|| exponent comes from a long range (in position space) interaction due to the 2pp 
singularity of the fermions. It drastically weakens the singularities associated with the transition. For example, in the 
weakly coupled spin liquid case the NMR T2 relaxation rate, which is given by Tj"^ ~ X]g(Xq)^j does not diverge as 
T — !■ 0, in contrast to transitions in 2D insulating magnets, where Tj"^ ~ 1/T. The T^^ rate is weakly divergent, see 



Eq. (3.16). The uniform susceptibility is nonvanishing as T ^ 0, but the leading temperature dependence is a power 
of T between and 1; see Eq. ( 3.19| ). 

In the strongly coupled spin liquid the scaling is again anisotropic. There is no singular dependence of x(Qi 0) on 
the direction of q. The only singular dependence is on the difference between q and "2pi^(0)", the vector spanning 
the Fermi surface and parallel to q. For cr > 1/3 the NMR relaxation rate {TiT)~^ is always divergent as T ^ (cf. 



Eq. (4.1)); diverges if cr > 1/2 (cf. Eq. (4.2)). Both quantities diverge more strongly as cr is increased. We also 



found that the leading temperature dependence of x is X '^"^ const + T^"*"^ where /3 is an independent exponent which 



tends to as cr becomes large and to 4/3 — 2cr as cr —> (cf. Eq. (3.19)). Our results for Ti, T2 and the uniform 
susceptibility are summarized in the Table. 

One reason for studying spin fiuctuations in two dimensional systems is that the high- Tc superconductors have been 
shown to have strong antiferromagnetic spin fluctuations. The two best studied high-Tc materials are La2-xSrxCu04^ 
and Y Ba2Cu^OT-s. In La2.xSrxCuOi, neutron scattering has observed peaks centered at incommensurate a:- 
dependent wavevectorso. At a; = 0.14 the incommensurability was shown to correspond to a "2pi^" vector of the 
LDA bandstructurda, suggesting that the results of the present paper should be relevant. NMR experiments have 
shown that the copper Ti rate ^^T^^ has the temperature dependence 



Cu 



(TiT)-i - 1/T (5.1) 
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for WOK <T < 500Jdla. The T2 rate has not-.been measured in this material, but in other high Tc materials with 
('^"TiT)"i - l/T, - T-'^ with i < x < lEZl. The uniform susceptibility is given by x - const + AT at least for 
150K < T < 4:00KtB. None of these properties are consistent with the weakly coupled spin liquid results of Section 
II. The strongly coupled spin liquid results with cr « 2/3 are in rough agreement with the data. 

In the YBa2Cu307-s materials, C'"(TiT)"i - 1/T and T^^ T""' for T greater than a 5-dependent "spin-gap" 
temperature and above this temperature x ^ const + A'T with A' ^-dependent. The neutron scattering indicates 
broad and flat-topped peaks centered at the commensurate wavevector (tt, tt). It is possible that the observed structuxa 
is due to several overlapping and unresolved singularities at a 2pF wavevectors. However, it has also been arguedllj 
that in Y Ba2Cu30j-s the magnetism is not driven by a 2pF instability, and is peaked at a commensurate wavevector 
because of a strongly peaked interaction. If the latter point of view is correct, the theory developed here is irrelevant. 
If the former point of view is correct, then we are forced to conclude that the relevant fixed point is the strongly 
coupled spin liquid with a « 2/3. 

One other aspect of the strongly coupled spin liquid requires further comment. We noted already that in this case 
X is singular at any Q spanning the Fermi surface. In a translationally invariant system this would predict peaks in 
the neutron scattering cross section on a ring of radius 2pF centered at the origin. As noted by Littlewood et alEl, 
for fermions on a lattice one obtains instead one or more curves traced out by the vectors Q connecting points with 
parallel tangents, and also one obtains additional families of curves displaced by reciprocal lattice vectors, G. One 
gets further peaking when members of different families of curves intersect. The resulting structure is very sensitive to 
the details of the Fermi surface and may resemble the data in Y Ba2Cu30T-s as well as the data on La2-xSrxCu04^. 
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TABLE I. Relaxation rates and susceptibility in different cases 





Weakly Coupled 
Spin Liquid (l/6<a<l/3) 


Strongly Coupled 
Spin Liquid (cr>l/3) 


i 

X{q ^ 0) 


I 

non-divergent 


1 

7'2<T-l/3 

1 

J'-JCT—I 

C + T'+f^ 
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